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ANALYTICAL METHODS 
FOR COMPUTING THE POLAR CURVES OF AIRPLANES. 



NOTATION. - The polar curve is in Cartesian coordinates, the 
curve representing the unit thrust y in terms of the unit 
drag x. 

If we express y and x in a parametrical form in terms 
of the incidence a , we note that within the usual limits of 
utilization, y is practically a linear function and x a 
parabolical function of a . 

The following computations will he made with this approx- 
imation. 

In the "Aerophile" of May 15, 1920, Mr. Rateau assumed, 
for the complete airplane equations in the form: 

( Y - Y 0 (l+ t) a ) 
( X = X Q (1+ I a 2 ) 
and defined a "typical airplane" for which 

Y 0 = 0.5897 n = 0.3 

X 0 = 0.1204 I = 0.017 

It seems to us interesting to give a method by which 
these equations can be established rapidly, taking as a basis 
profile tests made in the laboratories. 

We must, however, assume units differing from those of 
Mr. Rateau. 

1st. As units of Lift and' Drift we take y and x, the 
coefficients which the Germans call respectively C a and C w , 
which are 1,600 times greater than the Eiffel Ky and K x . 
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For that, it is sufficient to correct, for each value of 
lift, y, the measured values of x and a. 

These corrections represent induced drag and deviation, 
and ar^ given by Munk 1 s formulas. 

The profile polar curve is then considered as identical 
tc a parabola within the usual limits of lift, SO to 120, 
(corresponding to the Eiffel -Ky values varying from.. 0.0125 
to 0.075) . 

(3) x = x 0 + \ (y - y 0 ) 2 

In the same way, lift in terms of the incidence, is compared 
to a straight line: 

(4) a = p '+ p y 

P being the angle of zero lift, 

The G&ttingen laboratory has made the computations for the 
above 5 coefficients for 200 different profiles and has dram 
up a summary in the form of a Table, completed by two synopti- 
cal graphs. (Technische Berichte, II, 3; p,456.) 

A3 an example, we give a few profiles from this Table, cor 
responding to known planes, (See Table I, next page, 

W is the camber of the dorsal profile, D the maximum 
thickness, and H the thickness at 2/3 of the chord. 

These three numbers are given in fo of the chord. 



w If » 
TABLE I 



: : Characteristics x=x 0 fr\(y-y 0 ) : °--fe+P y.: 

No. of: :::::: : 

Prof.: Specification: : : : : : : 

: : W : D : H : y Q : x Q :10000\: p : 100P 



100 


: Sop with 


7.8: 


6.7 : 


5.4:45.3:1.50: 


4.5 


-3.3 


: 9.3 


135 


:Hanspeter : 


13.1:13.0: 


5.8:78.8 :1.80: 


3.3 


: -8,3 


: 9.8 


136 


:Gerhardt : 


7.3: 


7.1: 


3.4:46.8:1.36:11.3 


: -4.7 


:14.4 






16.3:16.3:11.8:50.6:3.16: 


5.4 


: -7.0 


:10.3 


357 


:Ago C IV 


: 9.5: 


8.5: 


6.0:57.7:3.03: 


5.0 


: -4.7 


:10.0 


358 


•.Hydro H Bbg 


:10.3: 


5.8: 


3.5:88.3:3-54: 


6.5 


: -5.6 


r 9.7 


359 


: AEG CIV 


*. 8 • 3 i 


6.7: 


3.3:48.9:1.95: 


5.3 


: -3.8 


:10.5 


360 


:Voisin 


: 9.3: 


4.1: 


3.9:91.3:3.67: 


6.1 


: -6.0 


: 8.7 


364 


:Fr'hafen 53 RI 


:10.0: 


6.5: 


4.4:73.7:3.38: 


6.3 


: -6.3 


:10.8 


365 


: « 41 RI 


:11.4: 


7.4: 


3.5:89.4:3.86: 


6.5 


: -6.4 


: 9.7 


368 


: Rumple r CIV 
: (upper wing) 


i 8*8* 


• 

6.0: 


• • • 

3.6:60.3:3.07: 


5.6 


• 

: -4.9 


• 

•10.7 


369 


: Rumple r CIV 
: (lower wing) 


• ♦ 

: 9. si 


6.6 : . 


• * • 

• • * 

4.0:70.4:3.08: 


7.5 


: -5.8 


•11.3 


398 


:Foiker Dr 1 . 


:13.7: 


13.6: 


7.5:73.8:3.33: 


4.4 


: -6.9 


:10.3 


306 


:MVA Spad 


: 6.8: 


.5.4: 


4.4:48.1:1.35: 


7.3 


: -3.1 


: 8.7 
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POLAR CURVE OF ANY MONOPLANE WING. 

Induced drag, which is added to that of the profile, varies 
as the square of the lift. 

In the same way, -induced deviation is proportional to lift. 
By the formulas of Bets we may assume: 



(5) 



( A x = ^y* 
( 



( A a = p.' y 

p. and jj.' are given by Bets in terms of the area S, of the 
span b and of tt . 



(6) 



H = 



1 
100 



1 — 



U-' = 



( 



(Kb)" 
S 



100 (K'bf 



K and K' are coefficients % 1 which depend on the geometri 
cal form of the wing. 

For rectangular wings K = K 1 = 1. 

In these conditions, the polar curve of the monoplane wing 
becomes 



(7) 



( 



x = ( \ + u. ) (y - ~-^—)+ x 0 + 

( X + \x ) \ + |i 



y 0 c 



a = ( p + ) y + P 

EXAMPLE. - For a rectangular wing of Aspect Ratio 
equations (6) give 



i 

6 



\x = 0.00053 



p.'- = 0.03 



For profile No. 100 (Sopwith wing): 

( x = 0.00098 (y - 20. 8) 2 + 1.998 

(8) ( 

( a = 0.123 y - 3.2 

It is interesting to compare the figures given by these 
formulas with those of the laboratories. 

1st. COMPARISON WITH GOTTINGEN. 



a 


:y computed :y measured :x 


computed: 


x measured 


0 


26 : 


25.3 


2.0 : 


2.1 


1.5 


38.2 : 


, 39.5 


2.3 


: 2.3 


3 


: 56.4 : 


52.2 : 


3.0 


: 3.1 


4.5 


: 62 . 6 : 


62.7 : 


3.7 


: 3.9 


6 


: 74.8 : 


74.7 : 


4.8 


4.8 


9 


: 99.2 : 


93.6 : 


7.2 


: 7.3 



We see that the approximation obtained is of the order 
the rfieasurements. 

The parabolic formula is thus justified. 

*'2niL COMPARISON WITH THE EIFFEL LABORATORY. 

Wing R A 40 of the Eiffel Laboratory is a Sopwith with 
an aspect ratio of 6. 

The test evaluated in coordinates x and y, gave: 

0° 5° 10° 

y • 23.8 58.6 87.4 

x 1.86 4.02 8.42 



a 8 ii 

For these same values of y, the values of x computed 
by formula (8) would be respectively 

1.99 3.40 6.35 • 

The test -results of the two laboratories are therefore 
not in agreement. 

The aspect ratio being the same, we think that these 
divergences may be attributed to a difference of thickness in 
the profiles. 

We have, in fact, for W, D, and H I 
Eiffel Profile: 7.3 - 6.4 - 5.5 

Gottingen Profile: 7.8 - 6.7 - 5.5 

The German profile being thicker, its fineness ratio is 
less for small angles and greater for large ones. 

We had the idea of computing the characteristic coeffic- 
ients of the polar curve of the Eiffel profile, and found the 
following values: 

y 0 = 3 B = - 3.8 

x Q = 1.51 P = 0.135 

\ = 0.00? 

The equation of the Eiffel model becomes: 

( x= 0.00143 (y - 19. 5) 2 + 1,83 

(9) ( 

( a = 0.155 y - 3.8 

This formula is the one we shall adopt in the following 
examples. 

i 



THE POLAR CURVE OF A CELLULE. 



There are several German methods for computing the induced 
resistances of a cellule based on the profile resistance of a 
monoplane wing. 

We will apply these methods one after the other to a eel- 
lule tested at the Eiffel Laboratory: 

Cellule RB9, Sop with profile composed of two rectangular 
and equal wings measuring 900 x 150. 
Gap = 160. Stagger (rear) =40. 

A) MUNK'S METHOD. - This method generalizes the formulas 
of Betz (5) and (6) . 



grams in terms: 

1st. Of the ratio £ of span to gap. . 

n 

2nd. Of the ratio S. of depth of stagger to depth of wing. 

t 

3rd. Of the angle e of difference between the angles of 




\i and p,' have the same form. 



The coefficients K K 1 and |i" are given by Hunk's dia- 



at t aclc • 



Here : 



bj. 
h 



900 
160 



5.6 



e 
t 



-40 



-0.266 



150 



c = 



K 



1.1 



1.05 
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( n = 0.0083 

( u' = 0.055 

( 

( u" = 0.4 • 
When all the computations are made, the equation of the- 
polar curve becomes: 

( x = 0.00176 (y - 15. 8) 2 + 1.93 

(11) ( 

( a = 0.18 y - 3.4 

B) PRANDTL'S METHOD. - In this method no account is taken of 
stagger. 

In terms of £ = 0.178 
we have, in the case of two equal wings, the coefficient of 
induction 

a = e = 0.515 

1 + 5.3 § 

which corresponds to a K of Munk equal to 



K = /— = 1.15 
J 1 + °- 

(instead of 1.11 as before). 

We have |a = 0.0008 

The equation of the polar curve becomes: 
(12) X= 0.0017 (y - 16. 4f + 1.92. 

C) FORMULAS OF BLASIUS AND HAMBURGER. - These are more pre- 
cise in the sense that they make for each angle the correction 
of Drift and the correction of Lift due to interaction. 

In the particular case we are concerned with, the two 
planes have the same dimensions and the same profile. 



The formulas are considerably simplified and give: 



(13) 



A y = _L_ ( cp ao + cpau) y 
3 h 



A x = — — ( cpwo + cpwu) x 
3 h 



The 4 functions cp are given by Tables and graphs in 
terms of the incidence and of the angle of stagger p . 
("Technische Berichte, II, 3, p. 341) . 



We have 



a. 



150 



3 h 3 x 180 



= 0.47 



cp ao 
cpau 
^ wo 



-0.13 
-0.09 
+ 0.013 
-0.007 



-0.11 
-0.14 
+ 0.018 
-0 .010 



10 

-0.S6 
-0.30 
+ 0.034 
-0.017 



The values of x and y for the monoplane wing are com 
puted by formula (9) for the Eiffel Laboratory profile. 

COMPARISON OF RESULTS. - In Table II we compare the figures 
given by each of our three computations with those measured 
by the Eiffel balance. 
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TABLE II. 











MUNX 


: PRANDTL 


BLASIUS 


:EIBFEL TEST 






( 


y : 


18 .9 


: 18.9 : 


31.2 


: 20 . 3 


a 


= 0 


i 
















( 


x : 


1.9 


1.9 : 


1.9 


1.9 






( 


y 


: 46.6 


: 46.6 


51.7 


: 49.1 


a 


- 5 


( 
















( 


X 


3.6 


: 3.5 


4.0 


: 3.9 






( 


y 


: 74.4 


74.4 


: 76.7 


: 74.4 


a 


= 10 


( 














( 


X 


: 7.6 


*: 7.6 


: 8.4 


: 8.1 



i7e see that these formulas are satisfactory for computing 
the parabolic polar curve of a rough draft, liunk 1 s method is 
rather to be recommended when there is a good deal of stagger 
and much difference between the angles of attack of the two 
planes of the cellule* Prandtl r s method is very suitable for 
straight cellules of unequal span. 

The method of Blasius is rather longer, but more exact, 
and takes into account the coefficients cp of quantities neg- 
lected by Munk and Prandtl. This method is more suitable for 
a thorough study of a chosen cellule. 

Finally, we see by the example of our Sopwith wing, that 
a slight difference of dimensions in the profile may cause much 
greater divergences than those occurring between the different 
methods. 

• POLAR CURVE OF AN AIRPLANE . 

To have the equation for the complete airplane, we have 
only to add together the drag of the cellule and the drag 
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t 

caused by structural resistances (landing chassis, rigging, 

bracing, empennages) . 

In theory, this drag should be separated into: 

1st. A profile drag which will vary, though very slightly, 
with the incidence. 

2nd. An induced drag due to the interaction of all these 
organs, with respect to each other and also with re- 
spect to the wing. 

Practically we may combine all these different kinds of 
drag in one constant coefficient. 

This approximation is all the more justified as what we 
are especially seeking to establish here is a typical equation 
on which we shall try successively all the combinations stud- 
ied in a rough draft. 

The tested values will, therefore, come in only in a rel- " 
ative fashion . 

We thus assume as parabolic equation of the glider: 

( x = ( \ + ^ )y 2 - 3 ^ y 0 y + x 0 + \ y| + v 
(14) ( 

( a = ( p + \x % )y + p + \x n 

MINIMUM DRAG. - This occurs for lift: 



and its value is 



(15) X n = x n + v + X ^ y 2 

The parabola of the airplane may be put in the form 

(16) x = ( X + |jl ) (y - Y 0 ) 2 + X Q 
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an equation of the same form as (3). 

For instance, if we utilize an airplane whose lift is 
just Y 0 , we see that maximum speed will be realized for X Q 
minimum. 

a) We reduce x 0 by taking a thinner .wing. 

b) We reduce v by diminishing structural resistances 
(fuselages, rigging, bracing wires). 

c) We reduce the last term, which may be written \x y 0 Y Q 
by diminishing induced resistances. (Action of the gap, as- 
pect ratio and stagger.) 

In reality, the requirements of the strength of materials 
necessitate a compromise which must be solved by tentative 
methods, and the problem is always brought back to a question 
of balance. 

FINENESS RATIO. - Maximum fineness ratio occurs for the 

lift 

/ x + X y 3 + v 

v = / o j_0 

V \ + \1 

and its value is 

(17) f = 1l = 2( .\ + n )Y r - 2 \ y Q 

On these formulas we can study direct the effect of the 
variation of one of the characteristic coefficients. 
MINIMUM POWER. - This is obtained by nullifying the derivate 
with respect to y in — £ — . It occurs for 



*r 15 *■ 



We have thus a simple relation between the three values 
of lift corresponding to minimum drag, maximum fineness-ratio, 
and minimum power. 

This equation is homogeneous in y. 

7 

We shall therefore have an analogous equation in a. 

In any airplane whatever, call: 

P -the angle of zero thrust. 

a Q the angle of minimum drag. 

a t the optimum angle, 

a z the angle of minimum power. 

These four angles verify the simple relation 

(19) ( a 3 - (3 f + 3( a 0 - p ) ( a 3 - p ) - 3 ( ^ - p f « 0 
independent of the origin chosen for the angles, and containing 
no coefficient of the airplane. 

EXAMPLE. .Take again Mr. Rateau 1 s typical airplane with its 
conventional sign: a. Q = 0 

Assume, as known: p = -3.33° and = 5.03°. 

The angle of minimum power, root of equation (19) gives 
us a =8.8° as computed by Mr. Rateau. 

2 

APPLICATION OF THE METHOD. 

Suppose we have an airplane having the cellule previously 
computed (equation 11) and structural resistance of 0.0016 kg, 
per m/sec. and per square meter. 

The coefficient v is 1.600 times greater. We have suc- 
cessively the following results: 
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CHARACTERISTIC COEFFICIENTS OF THE AIRPLANE. 



y 0 = 31 

x 0 = 1.51 
X = 0.0009 



p. = 0.00086 
\i' = 0.055 
v = 2.6 



6 = -3.8 
p = 0.125 
\l* = 0.4 



EQUATIONS OF THE POLAR CURVE 



x = 0.00176 ya - 0.0558 y + 4.97 
a = 0.18 y - 3.4 



MINIMUM DRAG. 



( Y n = 15.-85 
( 

( X 0 = 4.53 



a = -0.55' 
o 



The equation of the polar curve is written: 
x = 0.00176(y - 15.85)?- + 4.53 



MAXIMUM FINENESS RATIO. 



Y = 53.3 
f = 0.13 



a x = 6.18' 



MINIMUM POWER. 



Y 2 = 77.8 



a =10.6° 

2 



APPLICATION TO A CONCRETE CASE. - What engine will enable this 
airplane with 50 square meters of aerofoil to carry 1,500 kg. 
and fly at 173 km/hr. at an altitude of 3,000 m.? 
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In these conditions: Weight per square meter = 30. 
Static Pressure: q =105. 

We must have a lift: y = 100 j|| = 28.6 

The equation of the polar curve gives x = 4.8 
Total Drag: 

R = loo Sq - = ioo x 50 x 105 = 252 kg * 

Work per second: 353 x 43 = 13,300 kilogrameters. 
Assume a propeller efficiency 0.75 
The ratio of density = 0.74 
The nominal power required is 

T _ 12200 _ 300 horsepower. 

0 75 x 0.75 x 0.74 

CONCLUSION . 

This method is at least as precise as graphical methods. 
It is more rapid. 

Also, by this method all laboratory results can be carried 
in a pocket book in the form of Tables of Coefficients. 

Knowing the wind tunnel test of a wing and the performances 
of an airplane of the same profile, it is easy to verify the 
characteristic coefficients and, at the same time, the methods 
determining induced resistances. If the results had not been' 
previously checked, this method, entirely analytical, would 
enable us to discover the error. 

It is sufficient to vary the different data successively 
in close experimentation and to study their influence on each 
of the characteristic coefficients. 



Another important conclusion may be drawn from this art- 
icle. 

All French and foreign laboratories have already carried 
out many systematic tests- 

♦ 

But results cannot be compared, because comparative cali- 
bration has never been attempted. 

It would, however, be very simple to build two or three 
typical models of wings and cellules and to have such standard 
models in circulation amongst the various wind tunnels through- 
out the world. 

Only when that is done will experimenters everywhere be 
able to understand each other. 



